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1. Introduction 

The zero-range process (ZRP) is a simple lattice-based model in which particles move 
randomly to neighbouring lattice sites with rates depending only on the occupation of 
the departure site. This model was first introduced by Spitzer in 1970 pQ; more recent 
interest results from the observation that, for particular choices of hopping parameters 
the model exhibits a condensation phenomenon [21 El lU 03 H3 where a macroscopic 
proportion of particles accumulate on a single site (a real-space analogue of Bose- 
Einstein condensation). Condensation transitions occur in a wide variety of physical 
and non-physical contexts [3 El El EH HH U21 UH| and analytical studies of "toy" models, 
such as the ZRP, can provide important insight. The one-dimensional (ID) zero-range 
process has additional significance since it can be mapped to well-studied exclusion 
processes [Til ITo] . Condensation in the ZRP then corresponds to phase separation in 
the exclusion process [16J. For detailed discussion of recent progress on the ZRP and 
related models see, e.g., the review by Evans and Hanney [T7j . 

Most of these previous works study the ZRP defined on a ID periodic or infinite 
lattice. In a recent paper Levine et al. pH] considered the case with open boundary 
conditions, deriving the steady-state distribution and characterizing the condensation 
which occurs for strong boundary drive. Our contribution is to study the large-time limit 
of the current fluctuations through such an open ID ZRP with fixed initial configuration. 
This is in the spirit of previous studies of current fluctuations in the exclusion process ^H] 
and more general boundary-driven lattice gases |2()j . 

Specifically, the aim is to give the integrated current distribution and hence derive 
the large deviation properties of the measured time-averaged current. This helps to give 
insight into the nature of non-equilibrium steady states where large deviation functionals 
play a role similar to the free energy of equilibrium systems [21] . For example, there has 
been much recent interest in "fluctuation theorems" [2*2*1 EB1 EH EH1 EE] • By considering 
the dynamics under time-reversal, such theorems characterize the possible fluctuations 
of driven systems beyond the linear response regime and impose a symmetry property 
(to be discussed below) on the large deviation function. Essentially, this symmetry 
relates the probability of observing a given rate of entropy increase to the probability 
of observing the same rate of entropy decrease, leading to verifiable predictions for 
simulation |2Z] and experiment |28j . 

For the zero-range model considered here, our main physical finding is that, even for 
parameter values which result in a well-defined stationary state (i.e., no condensation), 
the current fluctuations undergo a change in character above some critical current. We 
interpret this qualitatively as a result of the temporary build-up of particles on some 
site(s) and argue that it is a precursor of the condensation transition. For unidirectional 
particle-hopping we show that the current distributions become spatially inhomogeneous 
for high currents. In this totally asymmetric case we also explicitly demonstrate the link 
to current fluctuations in the totally asymmetric simple exclusion process which we solve 
via the Bethe ansatz with particle-dependent hopping rates |29j. 
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Figure 1. Schematic representation of the ZRP on an open L-site lattice 

The remainder of the paper is structured as follows. In section |21 we define the zero- 
range model and summarize known steady-state results [IB] within the framework of the 
quantum Hamiltonian formalism. Section |3] contains a calculation of the integrated 
current distribution for the general ZRP (arbitrary hopping rates and boundary 
parameters) together with a discussion of the regime of its validity. In section 0] we gain 
deeper understanding by considering in detail the behaviour of current fluctuations in 
the totally asymmetric zero-range process. Then, in section EJ we present a determinant 
solution for the current distribution in this totally asymmetric model, obtained by 
utilizing Bethe ansatz results for the equivalent exclusion process. Finally, in section [H] 
we summarize our work and discuss some related issues. 



2. Zero-range process with open boundaries 

2.1. Definition of model 

We consider the ZRP defined on an ID open lattice of L sites — illustrated in figure ^ 
Each site I contains an arbitrary integer number of particles n which hop randomly 
to neighbouring sites (with exponentially distributed waiting time). In the bulk the 
topmost particle from site / moves to the right (left) with rate pw n (qw n ) where w n 
depends only on the occupation number n of the departure site and, by definition, 
wq = 0. Particles are injected onto site 1 (L) with rate a (5) and removed with rate 
•yw n (/3w n ). To facilitate later discussion of current fluctuations across different bonds 
we will label each bond by the site at its lefthand end, i.e., the Zth bond is between sites 
I and I + 1. 

For particular choices of w n , the model with periodic boundary conditions exhibits a 
condensation transition (see, e.g., [U 0111)1 HUB ED E2]) where, above some critical density 
p c , excess particles all accumulate on one randomly chosen site while the remaining 
sites have average density p c . In [18] it is demonstrated that (even for choices of w n 
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which do not give condensation on a ring geometry) open boundary conditions with 
strong driving can result in condensation at one or both of the boundary sites. In the 
present work, we first discuss current fluctuations for open boundary conditions with 
arbitrary w n (section |3J) before specializing to the totally asymmetric case with w n = 1 
(section 0]). This latter model can be mapped to the totally asymmetric exclusion 
process as discussed in section 03 

The zero-range process can be conveniently represented using the quantum 
Hamiltonian formalism [HH]. In this approach one works in the vector space (C 00 ) ^, 
defining a probability vector \P) = ^ n -^™l n ) with \n) the basis vector associated with 
the particle configuration n = (ni,ri2, . . . ,ul) and P n the probability measure on the 
set of all such configurations. \P) obeys the normalization condition (s\P) = 1 where 



(s\ = J2 n ( n \ an d {n\n'} 
represented by the Master equation 

d_ 

dt 



Within this formalism the ZRP time evolution is 



\P(t)) = -H\P(t)) 



(1) 



with 
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(3) 



and d is a diagonal matrix with the (i, j)th element given by WiSij. 



2.2. Stationary state 

The steady state of a stochastic process corresponds to the ground state of the quantum 
Hamiltonian. The left eigenvector with zero eigenvalue is just (s\, i.e., the constant row 
vector (1, 1,1,... 1). The corresponding right eigenvector obeys 

H\P*) = 0. (4) 

In ^H] h is shown that, for the ZRP, the stationary distribution is given by a product 
measure 

\P*) = \p*) <g> |P*) g, . . . ® \p*) (5) 
where |P*) is the probability vector with components 



p*{m = n) = 
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with the empty product n = defined equal to 1. Here Zi is the local analogue of the 
grand-canonical partition function 

oo n 

Z l = Z(z l ) = Y,z?l[w- 1 (7) 

n=Q i=l 

and the fugacities zi are uniquely defined by the following stationarity condition 

pzi - qzi+i = a - 721 = (3z L - 8 = c (8) 
with solution [TH] 

_ [(« + *) (p - g) - a/9 + 7*] (f ) - 7^ + a/3 (§) ^ 

7 (p - q - p) + P{p - q + 7) (j) 
and mean steady-state current 

c=(p-g) ^ — (10) 

7 (p - q - 0) + P{p - q + 7) (|j 

These results will be important for the calculations of the following sections. Note that 
the existence of the steady state, for given w n , is determined by the radius of convergence 
of Z. For parameters resulting in some Z\ outside this radius of convergence, a growing 
condensate occurs 



3. Current fluctuations 

3.1. Definitions and notation 

We are interested in the probability distribution of the total integrated current J/(i), 
i.e., the net number of particle jumps between sites I and I + 1 in the time interval 
[0,t]. It turns out to be convenient to study the generating function (e _AJi( ^) where 
the angled brackets denote an average over histories. In particular, we wish to find an 
explicit expression in terms of system parameters for the quantity e^(A) defined by 

ej(A) = lim -iln(e- AJiW ). (11) 

This characterizes the asymptotic current distribution which, for an ergodic system, is 
not expected to depend on the choice of initial particle configuration. 

As observed in |24j . equation (fTT|) implies a large deviation property for the 
probability distribution, Pi(j,t) = Prob(ji = j, t), of the observed "average" current 
ji = Ji/t. The long-time limiting behaviour is given by 

p,0\f)~e- | *&) (12) 

where ei{j) is the Legendre transformation of ej(A) § 

e t (j) = max{ej(A) - Xj}. (13) 

A 

§ This can be shown by performing a Fourier transform on the generating function and then using a 
saddle-point approximation to evaluate the resulting integral for large times. 
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We use the notation e/(A), ei(j) and Pi(j,t) to indicate the spatial dependence of 
the current distribution. It is clear that for a well-defined steady state the mean current 
must be spatially constant, i.e, (ji) = c for all I, but it does not necessarily follow that 
the full distribution of asymptotic current fluctuations is translationally invariant. We 
will return to this point later. 



3.2. Integrated distribution of input current 

Let us start by considering current fluctuations into the system, i.e., across the bond 
between the "reservoir site" and the first proper site 1. We represent the total 
integrated current across the "Oth" bond by the diagonal matrix Jq (whose diagonal 
elements are the set of all integer values) which acts on a state with basis vectors 
corresponding to the actual value of the integrated current across this bond. Then we 
introduce raising and lowering operators X ± on this space. These obey the commutation 
relation [X ± , Jo] = ZS ^X ± which implies 

e -\Jo X ± e xJ = e ^ x ± . (14) 

In the joint configuration/current state space the Hamiltonian is given by 

L-1 

H o = -{Yl H a 7 a ti - d i) + q( a i~ a 7+i ~ 4h)] 
i=i 

+a(X + a+ - 1) + 7(X~a^ - d x ) + 5(a+ - 1) + f3(a L - d L )|. (15) 

Now we wish to measure the expectation value 

(e- AJoW ) = (a, s|e- AJo e- Ho *|P, 0) (16) 

where the extra state labels refer to the current space and \P) is an arbitrary specific 
initial particle configuration obeying the normalization condition (s\P) = 1. Using the 
transformation (fT4*|) one obtains 

(e- AJoW ) = (s|e-^°*|P) (17) 



with 



L-1 



H = -{Y1 b( a * a m ~ dl "> + ^ a t a i+i ~ 



i=i 



^(a+e" A - 1) + 7(% e A - d x ) + 5(at - 1) + /?(a£ - d L )}. (18) 

In other words, the terms in the original Master equation which give a unit 
increase/decrease in the integrated input current are simply multiplied by e TA . We 
are particularly concerned with the long-time limiting behaviour which is given by 

lim(e- AJoW ) = (s|0)(0|P)e- go(A) * (19) 

t— >oo 

where eo(A) is the lowest eigenvalue of the modified non-stochastic Hamiltonian (fT%|) and 
|0) is the corresponding ground-state eigenvector. If the time- independent prefactor 
on the righthand side of equation (|TH|) is finite (meaning that the ground state is 
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normalizable in the sense that (0 1 0) and (s|0) are both finite) then we can identify 
e (A) with e (A), as defined by equation (fTTj) . 

To determine eo(A) for this case we make the ansatz that the ground state is still a 
product state of the form ©-(JTJ) but with modified fugacities and then use the following 
properties of creation and annihilation operators 

a+\0) = z^O) (20) 
aT|0) = *|0) (21) 
which can be explicitly checked with the matrices given in (jSJ). Then |0) satisfies 

L-l 

- H \0) = \^2(pz t - qz l+1 ){zf+ 1 di +1 - z^di) 
i=i 

+(ae~ x - -fz 1 )z^ 1 d 1 + (5 - f3z L )zl l d L - a + je x Zl - 5 + (3z L |0). (22) 

It is clear that if |0) is an eigenstate the coefficients of the gPs must cancel, leading to 
the condition 

pzi - qzi + i = ae~ x - ~jZ\ = /3z L - 5. (23) 

Comparison with equation (JHJ) shows that the required eigenstate of H is just given 
by the stationary product state of the ZRP with the replacement a — > ae~ x , i.e., the 
fugacities are 

[(ae~ A + 5){p-q)- ape~ x + j6) f 2 V _1 -76 + a(3e~ x (A 1 

zi = — —1 — • ( 24 ) 

l(p-q-/3)+(3(p-q + 7) 

The eigenvalue of this state is given by 

e (A) = -a + 7e A Zi - 5 + (3z L (25) 

and using (J2"2j) then yields 



(p-g)(e A - 1) 

e (A) 



a/3 (f) L_1 e- A -7^ 



7 (p - q - (3) + P(p - q + 7) (|) 



L-l 



(26) 



Since Hq is not stochastic the left ground state (0| is not necessarily (s\. In fact, 
using analogous methods to those employed above we can show that, for any choice of 
w n , it is a product state with (0j| = (l,zi,zf,...) where 

/3 7 (e A - 1) (| ) L ~ l + ie x (p - q - (3) + (3(p - q + 7) (j)^' 

7 (p - q - (3) + (3(p - q + 7) (j) 

As A — > 0, we see that |0) — > \P*), (0| — > (s\ and eo(A) — > 0, as expected. To 
argue that this product measure is indeed the required lowest eigenstate of H for all 
A we appeal to Perron- Frobenius theory ESI- We note that H contains positive 
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diagonal elements together with non-positive off-diagonal elements. It follows that the 
time evolution operator e~ Hot is a non-negative irreducible matrix which must have a 
real maximal eigenvalue corresponding to its only non-negative eigenvector. Now, it is 
straightforward to show that for all parameter values z\ > (and z\ > 0), so we have 
here found the positive eigenvector and (J2H|) must therefore be the lowest eigenvalue of 
H . However, note that this argument is not mathematically rigorous for an infinite- 
dimensional matrix. 

Equation for eo(A) is the chief result of this section and after Legendre 

transformation we obtain the rather cumbersome explicit expression for eo(j): 

{p-q)[aP{p/q) L - 1 + j5} 



lip-Q- P) + P(p-q + i)ip/q) 




L-l 



Aa(3 1 5{p/q) L -\p - q) 2 



hip - q - P) + P{v - q + i){p/q) L x ] 2 

2af3(p/q) L - 1 (p-q) 



q-(3)+/3(j)-q + 7)(p/q) 



L-l 



r + 



Aa(3^5{p/q) L - l {p - q) 2 



[21 



[l{p ~q-(3)+(3(p-q + l){p/q) L - l \ 2 _ 

In the following subsection we further discuss the regime of validity of (|2*Bj) , but first we 
make some comments on the form of our result. 

For the bulk symmetric case with driving boundary conditions taking the limit 
p = q = 1 in (J26|) gives 



(e* - l)(e-* tt /3 - 7 S) 

e » (A) = , + J + ' (29) 

Van Wijland and Racz [20] previously calculated the leading order term in the L — > oo 
expansion of this result by field theoretic methods. It has also been obtained by 
Bodineau and Derrida [36J via an additivity principle. 

As expected, our expression for e (A) obeys left-right symmetry (i.e., it is 
invariant under the transformation a <-> 5, /3 <->• 7, p ^ q and A <-> —A). Furthermore, 
it also demonstrates the Gallavotti-Cohen symmetry property. For a proof of the 
application of this fluctuation theorem [221EZI to stochastic particle systems see [2HEH]. 
In the present context it arises from a particular symmetry of the Hamiltonian Hq which 
we now outline. 

First, note from (fTU|) that the detailed balance condition for equilibrium (no steady- 
state current) is 

3 (if- 1 - (30) 

This corresponds to an equilibrium steady state \P* ) with fugacity z cq j = (5/ P)(p/q) l ~ L - 
In the non-equilibrium situation (|3Uj) is no longer satisfied but we can always write 
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a = a 



given by 



e E and 7 = 7'e E where detailed balance is obeyed for a', 7', (3, 5, p, q and E is 



e 2E = ^ ( P -\ • (31) 

Physically, one thinks of E as an effective external field imposed on the equilibrium 
model H2J. 

With these definitions we can show that H of equation (|18jl has the following 
symmetry property, 

H?(2E - A) = (P* q y 1 H (X)P* q (32) 

where P* is the diagonal matrix with equilibrium probabilities on the diagonal. 
Equation (J52*Jl means that the eigenvalues of Ho(\) and H Q (2E — A) are identical 
and hence, if the ground state of Hq(X) is normalizable, imposes the Gallavotti-Cohen 
relation 



L-l 



e (A) = e (2E - A) (33) 

which can be straightforwardly verified on the expression (|2fij). The oft-quoted form of 
the fluctuation theorem, 

PltlA „ e ~2Ejt (34) 

then follows via equations (fT2"j) and ((THJ). Physically, (p^jl relates the probability of 
observing a backward current with the probability of observing a current of the same 
magnitude in the forward direction. Unfortunately, such a relationship is difficult to 
check by simulation (or experiment) because of the exponentially small probabilities. 
Nevertheless, there has been some success in this regard for other models |2Z1 ESI EH] • As 
an aside, we note that the Gallavotti-Cohen fluctuation theorem reduces at equilibrium 
to the Green-Kubo formula and Onsager reciprocity relations — see jSHl 123 for details. 

3.3. Validity and spatial dependence 

The expressions (|2%j) -(|27 j) are general results independent of the hopping rates w n . 
However, the form of w n enters in the normalization of the ground state and hence the 
regime of validity of expression (|26J) for eo(A). The normalizing factor for (0 1 0) is given 
by the generalization of (J7J) 

00 n 

71=0 1=1 

with fugacities given by ()24)1 and (|27jh We note that, even for parameters for which 
the stationary state exists and there is no condensation, there may be only a restricted 
range of A for which (|33|) is convergent and (s\0) is also finite. The characterization and 
physical description of current fluctuations outside this regime is a central goal of this 
paper. 
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A related issue is the important question of whether the current fluctuations are 
the same across all bonds of the system, i.e., whether e;(A) is independent of I. Let us 
consider measuring the current between sites / and l + l in the bulk of the system. Then 

( e -M(*)) = ( 8 \e- 6lt \P) (36) 

with 

L-l 

Hi = -{ Yl [p( a k a t+i - d k) + <i( a i a k+i - 

+ [p(e~ x a^af +1 - di) + g(e A a+a m - d i+1 ) 

+a(a+ - 1) + 7(«r - d i) + 5(a+ - 1) + p{a~ L - d L ) j (37) 
Now, we define the operator 

Y r i = e AE *=i d *, (38) 
and note the generic integrated commutation relation 

for X(k) an arbitrary function of k. Using this commutation property we can show 

( e -M(*)} = (slY^-^Yr^P) (40) 
and hence we find, with the lowest eigenvalue eo(A) of H , the long-time limit 

lim(e- AJi(t) ) = (s|Y,|0)(0|YT 1 | J P)e" eo(A) *. (41) 



t— >oo 



Here, as before, |0) is the ground-state eigenvector of H and if the prefactor is finite 
(which now requires the additional condition that (s|YJ|0) is finite) then the current 
fluctuations across all bonds are seen to be the same (to leading order) as those going 
into the system, i.e., ej(A) = e (A) for all I. We have verified this by explicit calculation 
of ei(A) and ez,(A). 

Thus we conclude that our result is guaranteed to describe current fluctuations 
across any bond I but only for the range of A for which (s|0), (0 1 0) and (s|Y/|0) are all 
finite. From the form of the fugacities ()24j1 . (|27j) and the normalizing factors, we see 
that the appropriate range of A will depend both on the rate parameters (i.e., a, (3, 7, 
S, p, q) and on the details of w n . In order to further understand the significance of this 
regime of validity, and the behaviour of current fluctuations outside it, we examine in 
the next section the specific example of the totally asymmetric zero range process with 
w„ = 1. 



4. Totally asymmetric zero-range process 

4-1. Model and motivation 

We here consider the totally asymmetric version of the model defined in section viz. 
7 = 5 = q = with p = 1 for normalization of time. The choice w n = 1 + b/n is a generic 
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a 




1 

Figure 2. Single-site TAZRP with w n = 1 



model which, for periodic boundary conditions, gives condensation at high particle 
density for b > 2 (U |5J [TBI USB EH E21- For open boundary conditions condensation 
can occur even for b < 2 |18j : here we study the case 6 = for simplicity. 
For this model, from equation @, the stationary fugacities are given by 

zi =a for l^L (42) 

n = j (43) 

For a normalizable steady state one must have a < 1 and (3 > a; other cases result in 
the occurrence of a boundary condensate [18^ . Our aim is to study current fluctuations 
for all X, j in the case where the steady state is well-defined. For definiteness we mainly 
consider a < (3 <\. 

We first study analytically the toy case of a single site (section 14.2(1 before 
presenting heuristic arguments for the behaviour of current fluctuations in larger systems 
(section l4.3p . These arguments will be later confirmed by the Bethe ansatz approach of 
section [SJ 

4-2. Current fluctuations in single-site case 

It is illuminating to study the single-site (two-bond) TAZRP for which explicit 
calculation of e(A) is possible for all A. The model is illustrated in figure we choose 
(3 > a to ensure a well-defined steady state. 

Our main results are that the current fluctuations into the system obey, 

Po(j,t) ~ e - tla - j+jln{j/a)] (44) 

for all currents, while those out change in character at j = (3, viz., 

{ e -t[a-j+jin(j/a)] j < (3 

e -t[a-j+jln(j/a)] x e -t\P-j+jWj/ffl j > ^ 

The remainder of the subsection is devoted to the derivation of these expressions. 
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Let us first recap the results of section |3] as applied to this case. Here H of 
equation (fTHj) is given by 

H = - [a(afe- x - 1) + fife - di)] (46) 
which has ground-state eigenvalue 

e (A) = a(l-e- A ) (47) 
leading to a current distribution with 



eo(j) 



a 



j + j In 



J 



a 



(4* 



This corresponds, as should be expected, to a Poissonian distribution of particle jumps 
from the reservoir onto the site. According to the argument of section E?"2*l this result is 
only guaranteed to be valid if (s\0) and (0|0) are both finite. However, since the current 
entering the system is independent of the site occupation, the Poissonian distribution 
must hold for all j and hence in this case eo(A) is given by (}4*Tj) for all A. 

As we saw in section EHfl the regime in which ei(A) must be equal to e (A) (in 
other words, current fluctuations across the outgoing bond have to be given by the 
same Poissonian distribution) is determined by the range of A for which (s\0), (0 1 0) and 
(s|Yj|0) are all finite. Now, for the totally asymmetric case we find from equation (J27|) 
that zi — 1 so that (0| = (s\ and the first two conditions are equivalent. Furthermore 
it turns out (since z% is here proportional to e~ A ) that this condition also guarantees 
that (s|Y/|0) is finite when the steady state exists. Hence, the range of validity of 
ei(A) = a(l — e _A ) is simply determined by the range of normalizablity of |0). 

In this single-site case, we find that |0) has fugacity Z\ = ae~ x /j3 so that, for w n = 1, 
the condition for normalizability is ae~ x < (3. Via the Legendre transformation of (JT3j) 
this gives a current condition j < (3. In order to determine the behaviour of current 
fluctuations out of the system when this condition is violated (i.e., ae~ A > /3, j > (3) 
we now explicitly calculate ei(A). In other words, we look for the lowest eigenvalue of 
the matrix 

H l = -[a(af - 1) + P(a^e~ x - d x )\ (49) 

for which the eigenvector is well-defined. || 

Using the representations of creation and annihilation operators (J3J) together with 
the diagonal matrix d\ we see that the right eigenvector corresponding to eigenvalue A 
satisfies 

/ a -0e- x . . . \ 



-f3e~ x 
a + (3 



a 

-a 





-pe~ x 
a + f3 
—a 






-/3e" A 
a + f3 
—a 



\ 



J 



( bo \ 

h 



( b \ 



A 



bo 

h 
bs 



(50) 



V ••• / 



|| The results of this subsection may also be obtained rigorously by calculating (s\e Hlt \P) in integral 
form and examining the long-time behaviour. 
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We thus obtain a mapping for the elements of the eigenvector 

-a6 n _i + (a + - A)b n 
b n+1 j^- x (51) 



with initial condition 

'a — A 

This mapping can be written as a transfer matrix equation 

b 



h = ( -ct ) b °- ( 52 ) 



(53) 



[ bn+1 1 










where the transfer matrix has eigenvalues 



t*± = • ( 54 ) 



a + (3 - A ± y/(a + (3 - A) 2 - 4a(3e~ 
2(3e~ 

Similarly the elements of the corresponding left eigenvector (a , a 1; a 2 , a 3 , . . .) of Hi 
are given by the mapping 

f a- A\ . . 

ai =1 — - — I a (55) 

-/3e" A a n _i + (a + /3 - A)a n 

a n+ i = (56) 

a 

corresponding to a transfer matrix with eigenvalues 

_ a + (3 - A ± y/(g + (3 - A) 2 - 4a/3e- A 

From the properties of these transfer matrix eigenvalues we can determine the large 
n behaviour of a n (and b n ). Specifically, the limit of the ratio a n /a n _i is given by the 
largest magnitude eigenvalue except in the case that the initial ratio ai/a coincides 
with the smaller eigenvalue. Now, if lim^oo |(a Tl 6 n )/(a n _i6 n _i)| > 1, then the series 
ao&o + a ibi + ^2^2 + • • • CL n b n • • • is divergent (D'Alembert ratio test) and hence the 
eigenvector of Hi with eigenvalue A cannot be normalized. Such an analysis leads 
to the conclusion that the lowest eigenvalue for which the corresponding eigenvector 
can be defined is given by^f 

fa(l-e _A ) ae~ x < (3 

ei(A) = < , (58) 

[a + (3- 2^/a(3e~ x ae~ x > (3. 

Using the Legendre transformation (|13|) this gives 

a. — j + j In — j < (3 

<dj-t)={ a j j (59) 

a - j + j In - + (3 - j + j In - j > (3. 

a (3 

% One can readily check that in both cases (s\l) is finite with |1) the right eigenvector corresponding 
to ei(A). 
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In other words, for currents greater than (3 the current fluctuations out of the site are 
different to those into the site. This behaviour can be explained heuristically as follows. 

For j < (3, the current into the site is smaller than the average attempt rate 
out. This means that hopping across the two bonds is not independent and there 
is no chance for particles to build-up on the site — all particles which arrive there 
move almost immediately out of the system. Hence, for j < f3, the limiting factor 
in observing a current fluctuation with j 1 = j is just the input current j , i.e, 
Prob(ji = j, t) ~ Prob(j > j, t) and we find that ei(j, t) = e (j, t) = a - j + j \n(j/a). 

In contrast, for j > (3, the current into the site is greater than the attempt 
rate out leading to a temporary build-up of particles which we dub an "instantaneous 
condensate". (We emphasize that for the parameters we study, a < (3, there is no 
permanent condensation, the temporary build-up of particles occurs only in those 
unlikely realizations where j > (3.) In this case, at large times the probability of 
finding at least one particle on the site is essentially unity, so the two bonds behave 
as independent Poissonian processes (with means a and (3 respectively). Hence, for 
j > (3, Prob(ji = j,t) is given by the product of Prob(j > j,t) with the independent 
Poissonian probability for observing a current j across the outgoing bond. This gives 
eiO', t)=a-j+j Hj/a) + [3 - j + j Hj/(3). 

4-3. Current fluctuations in L- site case 

Based on this picture of "instantaneous condensates" we propose that current 
fluctuations across the different bonds in an L-site TAZRP, with w n — 1 and a < (3 < 1, 
are given by the following expressions: 

• Input bond 

Po(j,t) ^e-^ a - j+j ^ j ^. (60) 

• Bulk bonds, 0, L 

{ e -t[a-j+j\n(j/a)] j < 1 

e - t [a-j+jln(j/ a )] x e - t (l-j+j In j)l ■ > ^ ( 61 ) 

Note the /-dependence here! 

• Output bond 

{ e -t[a-j+j\n(j/a)] j < (3 

e -tia- j+j HjM] x e -t[V>-j+3W/f>)] (3 < j < 1 (62) 

e -t[a-j+jln(j/a)} x Q -t{l-j+j In j)(L-l) x Q -t\J3-j+j Hl/P)} j > I. 

In the next section we argue that these expressions can be rigorously obtained via a 
Bethe ansatz approach. First we present here some general arguments for their form 
and significance. 
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The low current behaviour can be directly obtained from the analysis of section El 
The lowest eigenvalue of 

"L-l 

^^p(a^a i t fl — di) + a(afe~ x — 1) + /3(a^ — di) 



Hn = - 



1=1 



(63) 



is 

e (A) = a(l-e- A ) (64) 
for any L, again giving a Poissonian distribution of current into the system, i.e., 

eo(j) = a-j + jhi-. (65) 
a 

The ground state |0) of H corresponding to eigenvalue eo(A) diverges on site L for 
ae~ x > (3 (j > (3) and on all sites for ae~ x > 1 (j > 1). For current fluctuations above 
these limits we argue that the divergence of this ground state indicates the temporary 
build-up of particles — "instantaneous condensates" — in this sense it is a precursor of 
the condensation which would occur for (3 < a or a > 1. This causes the spatial 
homogeneity of the current fluctuations to be broken. In summary, we expect to find 3 
different regimes: 

I. A > \n(a/f3), j < f3: Current fluctuations across all bonds are identical with 
ei(X) = e (A) = a(l - e~ A ) and q(j) = e (j) = a -j + j ki(j/a) for all I. 

II. ln(a://3) > A > In a, /3 < j < 1: In this case the eigenvector corresponding 
to eo(A) is not normalizable for site L. Physically, we argue that for j > j3 the 
current into site L is larger than the attempt rate out leading to a temporary 
build-up of particles on site L. The probability of observing a current fluctuation 
with ji = j is then given by the probability of observing — j multiplied by 
the independent Poissonian probability (mean (3) of particles at site L exiting the 
system with rate j. For all other sites the ground state is well-defined, behaving 
like a system of size L — 1 with effective exit parameter f3' = 1. Hence, for / ^ L, 
ei(j) = e (j). 

III. A < lna, j > 1: Here the eigenvector corresponding to eo(A) is ill-defined on 
all sites and hence we expect "instantaneous condensates" on all sites. Each bond 
then acts as an independent Poisson process. A particular feature of this totally 
asymmetric model is that the current through a given bond is dependent only on 
the bonds to the left. Hence, for j > 1, Prob(ji = j, t) is a product of / + 1 
Poissonian probabilities and the distribution of current fluctuations through every 
bond is different. 

For fixed (3 < 1, this behaviour can be summarized in the 2D phase diagrams of figure El 



For the alternative case a <1 < (3 one expects instantaneous condensates on sites 
1 to L — 1 for j > 1 and on all sites for j > (3. The probability distribution for current 
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Figure 3. Phase diagram of TAZRP (with w n = 1) for case (3 < 1 in (a) a-e A plane 
and (b) a-j plane. Phases I, II and III correspond to the different regimes of current 
fluctuations described in the text. No steady state exists to the right of the solid line. 



fluctuations out of the system will be accordingly altered to 

{Q-t[at-j+j lnO'/ct)] 
e -t[a-j+j\n(j/a)] x e _t(l-j+j In j){L-l) 
e -t[a-j+jln(j/a)] x e ~t(l-j+jln j)(L-l) x e -t\j3-j+j ln(j//3)] 

5. Application of Bethe ansatz results 

5.1. Mapping to exclusion process 

In this section we demonstrate how the application of new Bethe ansatz results [2H] for 
the totally asymmetric exclusion process with particle-dependent hopping rates provides 
an alternative analytical approach to rigorously derive the TAZRP results of the previous 
section. In particular, we wish to recover equations ([DD]) -([D2 ]) for the site-dependent 
probability of observing current fluctuations with a given j. 

The general ZRP can be mapped to an exclusion model (in which lattice sites cannot 
be multiply occupied) by representing each bond as a particle and the occupation of the 
site between a pair of bonds as the number of interparticle vacancies. For example, the 
particle configuration shown in the L-site ZRP of figure ^ (with each bond labelled by 
its lefthand site) is equivalent to the exclusion picture of figure 0] with L + l particles — 
note that with our numbering convention the Zth bond of the ZRP maps to the (I + l)th 
particle in the exclusion process. The hopping of a ZRP particle in a rightward (leftward) 
direction across a given bond, then corresponds to the equivalent exclusion particle 
moving one lattice site left (right). For the case w n = 1 (the pure chipping process [3D] ) 
the hopping attempt rates are independent of the particle separation leading to the 



J<1 

1 < j < P (66) 
j>P. 
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1 2 3 4 5 6 L-l L L + l 

Figure 4. Exclusion model corresponding to zero-range process of figure ^ with 
w n = 1. Bulk particles hop to vacant nearest- neighbour sites on the left (right) with 
rate p (q). Particles 1 and L + l have different hopping rates corresponding to the 
boundary processes of the ZRP. 



(asymmetric) simple exclusion process with parameters assigned as in figure EJ 

We emphasize that the ZRP on a finite open lattice maps to an exclusion process 
with a fixed number of particles on an infinite lattice where the endmost particles 
have different hopping rates from the others. An initially empty lattice in the ZRP 
corresponds to the particles in the exclusion process all occupying adjacent sites (no 
interparticle vacancies). 

In particular, the open-boundary totally asymmetric zero range process (TAZRP) 
of section H] maps to the totally asymmetric exclusion process (TASEP) on an infinite 
lattice. The pure TASEP (particle-independent hopping rates) has been extensively 
studied, for both periodic and open boundary conditions, and the steady state solved 
by a variety of means including recursive techniques [H], a steady-state operator algebra 
formulation [12] and application of the Bethe ansatz to the corresponding quantum-spin 
system Recently Rakos and Schutz [2H] have extended the Bethe approach for 

calculating current fluctuations to the general case of particle-dependent hopping rates. 
We now summarize the result of that calculation, reordering the numbering of particles 
in order to make the connection with the TAZRP more transparent. 

Consider an infinite one-dimensional lattice containing N particles with left hopping 
rates v\,v 2 , . . . ,vn (from left to right). If the particles are initially located at 
sites yi, y 2 , ■ ■ ■ , Vn then the probability P(x,y,t) of finding them at time t on sites 
xi,X2, ■ ■ ■ ,xn (with lattice coordinate increasing to the left and x\ > yi V/) can be 
written in the determinant form 

N 

P(x,y,t) = Y[(e- tv % Xi - yi ) x (67) 

i=l 

F^iixx-y^t) F lj2 (x 2 - yi,t) ... F 1}N (x N - y h t) 
F 2 ,i(x 1 -y 2 ,t) F 2>2 (x 2 -y 2 ,t) ... F 2:N (x N - y 2 , t) 

F N>1 (xx-y N ,t) F N>2 {x 2 -y N ,t) ... F N>N (x N - y N , t) 
with F-functions defined as 

N N 

F r>s (x,t) = — fe^z*- 1 J] (1-Viz) J] (l-v t z)- l dz (69) 

i=r+l i=s+l 

where the integral is to be taken along a circle of radius e around the origin of the 
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complex plane. These functions obey the following identities 



.''2 



£ v s x F r>s (x,t) = v s ^F r>s ^(x u t) - v s ^ +l F r ^{x 2 + l,t) (70) 

x=x\ 

v r+1 x F r>s (x, t) = v r+1 xl F r+liS (x u t) - v r+1 X2+1 F r+hs (x 2 + 1, t) (71) 



x=x\ 



with important special cases for x\ = x 2 = x 

F r>s (x,t) =F P|J ,_ 1 (x,t)-v a F r>a _i(x + l,t) (72) 
F r>s {x, t) = F r+ltS (x, t) - v r+1 F r+1>s (x + l,t). (73) 

In the next subsection we show how the result can be applied to obtain 
the current fluctuations in the TAZRP. Note that for brevity we will write F(x) as 
a shorthand for F(x,t). 

5.2. Determinant solution for current fluctuations in TAZRP 

As outlined above, we map the L site TAZRP to an (L + l)-particle TASEP where the 
particles hop with rates v\ = a, vl+i = (3 and i>/^i,£+i = p = 1. To study current 
fluctuations across the Zth bond of the TAZRP we need to know about the number 
of jumps made by the corresponding exclusion particle. Importantly, for this totally 
asymmetric model, this is independent of all the particles behind and thus we need only 
consider the first l + l particles. For simplicity we will start from the initial configuration 
of adjacent exclusion particles (equivalent to an initially empty lattice in the ZRP); as 
previously discussed, we expect the long-time decay of the current fluctuations to be 
independent of the specific initial particle configuration. 

Now we can use the Bethe ansatz solution to obtain the distribution of the distance 
travelled by the (/ + l)th particle 

Prob(xi > x + I, X2 > x + I — 1, . . . , X[ > x + 1, Xi + \ = 

yi = I, y2 = l - l, • • • yi = 1, yi+i = o, t) 



x\ 



xi — '1 



£ £ ... £ P({x 1 ,x 2 ...x l ,x},y,t) 



(74) 



Xl=X+l X2=X+l — 1 Xi=X+l 

l+l 



e tVi Vi X ) x 



i=l 



Fi, {x) 
F 2 ,o(x + l) 



f 2 ,i(x) 



F l+1 ,„(a; + F l+1>1 {x + I - 1) 
i+i 

i=l 



• •• Fi,,_i(ar - Z + 1) F 1)l+1 {x-l) 
••• F n „r{x - I + 2) F 2 , l+1 {x - I + 1) 

• •• Fj + u_i(a; + l) F l+ljl+1 (x) 



(75) 
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F i+i,o( x ) 
F l+lfi (x + 1) 



Fi + i fi {x - 1) 
Fi+i,o{x) 



... F l+1>0 (x - I + 1) 
... F l+1;0 (x - I + 2) 



F l+1>l+1 (x - 1) 



(76) 



Fi+i,o(x + l) F l+ i t0 (x + 1 - 1) ... F/+i j0 (x + l) F i+ i ii+ i(s) 

where the intermediate steps use determinant manipulation together with F-function 
identities given in the preceding subsection 

To obtain an explicit expression for the current distribution across the Ith TAZRP 
bond we set x = J = jt and drop the redundant subscript to yield 
i+i 



Pi 



(j^)=ri e " 



■t(vi-jln Vi) 



X 



i=l 



D (jt,t) D (jt-l,t) 

D Q (jt + l,t) D (jt,t) 

D Q (jt + l,t) D Q (jt + l-l,t) 
with elements 



D (jt-l + l,t) D l+1 (jt-l,t) 
D (jt-l + 2,t) D l+1 (jt-l + l,t) 



D {jt + l,t) 



Di+i(jt,t) 



(77) 



D.(x,t) = — 

Z7T1 



l+l 



ViZ) 



^dz. 



(78 



i=s+l 



For the TAZRP with initially empty lattice this expression is exact for all times; 
we expect the long-time limit to give the initial-configuration-independent expressions 
suggested previously (|HU|) - (|fi2|) . 

To understand how different regimes of current fluctuations emerge from this 
determinant solution, let us examine the case of a two-site TAZRP (equivalent to a 
three-particle TASEP). As in section HJ we take a < (3 < 1. 

We begin by considering the current fluctuations across the Oth bond of the TAZRP. 
In this case the determinant is just a single integral 

1 ''' ' " ' (79) 



-t(a—j In a) 



2vri 



3*/V t - 1 dz 



and a saddle-point approximation leads straightforwardly to the asymptotic behaviour 

Po(j,*)~e-'( Q -^ ln ^ (80) 

in agreement with our previous arguments (recall that the current into the system is 
always described by a Poissonian distribution). 

For current fluctuations across the 1st bond of the TAZRP (i.e, between sites 1 and 
2) we need a 2 x 2 determinant 

_L ^e'/V-^l - z)-\l - az)- l dz ± §e t / z z^ t - 2 dz 



Pi(j,t) 



-t(a—j In ct+1) 



1) 



J_ ^ e */V(l - z)- 1 ^ - az)- x dz ± § e l / z z jt ~ l dz 

Again we argue that the long-time behaviour can be extracted by saddle-point analysis 
of the individual elements. However, more care is now needed because of the poles in 
the integrands of the D functions. Specifically, we note that these poles are at z = 1 
and z = a -1 while the saddle-point is located at z = 
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For j > 1, both poles lie outside the saddle-point contour. In order to obtain 
the leading-order non-zero term in the determinant it is necessary to expand the t- 
independent factors in the integrands to second order about the saddle-point. This 
gives the long-time behaviour 

Pi(j, *) ~ e- t[a ~ j+jHj/a)] x e -*0--3+3toj) m ( 82 ) 

In contrast, for j < 1, in deforming the contour through the saddle-point we 
include the pole at z — 1. Subtracting the residue at this pole from the saddle-point 
approximation for each D Q integral yields the changed leading-order behaviour for the 
determinant 

pi(j,t) ~ e -* [Q - i+iln{j ' /a)1 . (83) 

For j < a the pole at z = a -1 is also inside the contour. However it turns out that 
this results in the addition of a sub-leading term and the long-time limit is unchanged. 
Hence, for the full range of j, we recover the current fluctuation distribution predicted 
in section HI [cf. equation (pTTj) with 1=1]. 

To determine current fluctuations across the 2nd bond of our two-site TAZRP (i.e., 
the current fluctuations from site 2 out of the system), we must add a third particle 
with hopping rate f3 to our TASEP picture and evaluate the 3x3 determinant given 
by ()77j) with 1 = 2. Here a similar saddle-point analysis with careful treatment of the 
poles at z = a -1 , z = 1 and z = leads to 

{ e -t[a-j+jln(j/a)] j < (3 

e -t[a-j+jHj/a)] x e - t[ p- j+jH j/f3)] (3<j<l (84) 

e -t[a-j+jln(j/a)} x ^-tfl-j+jhxj) x e ~t[f3-j+j ln(j / (3)} j >1 

in agreement with equation ()62|) for the case L = 2. 

The validity of this saddle-point approach is confirmed by numerical evaluation of 
the determinants. In figure El we show the result for the second bond for several different 
times. Convergence to the long-time limit of (J54*|) is clearly seen. Good agreement is 
also shown with simulation results for current fluctuations in the totally asymmetric 
zero-range process with empty initial lattice. + For other initial configurations we found 
slightly different curves for small times but, as expected, the same limiting distribution 
for long times. Unfortunately, since the chance of seeing a current fluctuation away from 
the average becomes exponentially more unlikely with time, obtaining large deviation 
simulation data at long times is very computationally intensive. 

In figure |H1 we show numerical results for the current fluctuations across the three 
different bonds in our toy system. The results are in full agreement with the analytical 
expressions and graphically illustrate the spatial inhomogeneity of the fluctuations for 
large currents. 

It is obvious that to study current fluctuations in a TAZRP with L sites (L + 1 
bonds) requires an L+l particle TASEP and hence the evaluation of an (L + l) x (L + l) 

+ For efficient simulation we utilized the exact mapping of the TAZRP to a last passage percolation 
problem |45| . 
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J 

Figure 5. Current fluctuations across the 2nd bond (out of the system) in a 2-site 
TAZRP with a = 0.3, p = 1 and /3 = 0.5. Graph shows e 2 (j) = limt-^l/t) lnp 2 (j> *)■ 
Solid symbols on dashed lines result from numerical evaluation of the determinant 
solution at times t = 100, 300 and 1000. Crosses are simulation results for t = 100. 
For increasing times, distributions converge on the long-time limit of equation 18411 
shown by the solid line. 



determinant. This becomes increasingly tedious for large L. However, it is clear that 
each bulk particle with hopping rate one gives a factor (1 — z)^ 1 in the Do integrals of 
the determinant and hence that a saddle-point expansion for large times with careful 
treatment of the poles (as above) will lead to equations (jEIJ) and (JH2J)- 

We note that the exclusion process picture gives further qualitative explanation 
for the behaviour of the current fluctuations in the different regimes. One sees that 
for j < 1 the bulk exclusion particles all pile up behind the leading particle, whereas 
for j > 1 they move independently In other words, since particles cannot overtake, 
a bulk particle (with unit hopping rate) always has the same average current as the 
first particle if the latter is moving with speed less than one. On the other hand, for 
a particle to move with current greater than unity requires all the particles in front of 
it to be moving independently with j > 1. Similarly, the last particle in the TASEP is 
located immediately behind the penultimate particle for j < (3 but moves independently 
for j > p. 
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Figure 6. Current fluctuations across the three different bonds of a 2-site TAZRP with 
a = 0.3, p = 1 and [3 — 0.5. Graph shows e;(j) = limt_ >00 (l/t) hip/(j, t). Solid symbols 
result from numerical evaluation of the determinant solution at time t = 1000. Lines 
are the long-time limits given by equations l|80|l and H82I) -(I84 |1 . As expected, curves 
for the different bonds split at j = (3 and j = 1. 



5.3. Conjecture for disordered TAZRP 

Noting that the determinant solution (|77|) is valid for any choice of vi motivates us 
to consider the behaviour of current fluctuations in the totally asymmetric zero-range 
process with bond-disordered rates. We consider an L site TAZRP (again with w n = 1) 
where a particle moves from site / to site I + 1 with rate r\. The rate of hopping from 
the lefthand reservoir site into the system is denoted by r . The pure case considered 
above corresponds to the choice r = a, Tl = (3 and n, r 2 , . . . , Tl-i = 1- This problem 
can of course be represented as an (L + l)-particle exclusion process with hopping rates 

vi = 

One can easily show that the fugacities for the stationary state of the disordered 
TAZRP model are given by 

zi = - (85) 

n 

zi = 1, (86) 
and thus the condition for a well-defined stationary state is r < r t for all I. In fact, the 
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general conjecture to be presented below also describes current fluctuations in the case 
where this condition is violated and condensation occurs on one or more sites. 

In our analysis of the current fluctuations we are interested in the product state with 
fugacites obtained from (|85|) by the substitution r — > r e~ A . If this is to be normalizable 
on site I then we clearly require roe~ A < r; corresponding to the current condition j < rj. 
Since, in the totally asymmetric model, the current across a given bond is only dependent 
on the occupation of sites to the left it is clear that an important limiting role is played 
by the smallest hopping rate leftward from the bond under consideration. Hence, for 
convenience, we define 

rt = minjrfc}. (87) 

k<i 

Consideration of the form of the determinant solution (|77j) then leads to the 
following generalization of our pure results 



r i - j + j In — j < rj 



•/ j \ ( 88 ) 

22[r k -j+j\n—\@(j-r k ) j > r* 

v k=o ^ Tk ' 

where Q(x) is the Heaviside step function. Note, in particular, that the limiting 

behaviour of the current fluctuations is independent of the ordering of hopping rates 

leftward from the bond in question. In fact it follows from (J77J) that, for an empty initial 

condition, this is also true for all finite times. 

In other words, for currents less than the minimum hopping rate r ; * the current 

fluctuations are always controlled by the "slowest bond"; for higher currents there are 

"instantaneous condensates" on all sites with exit hopping rate r\ < j. 



6. Discussion 



In this paper we have considered the large-time limit of the current fluctuations in 
the zero-range process with open boundaries and some fixed initial configuration. 
For arbitrary boundary parameters and hopping rates we obtained a general result 
[equation (J2HJ)] for the integrated current distribution and discussed the regime of 
its validity. We then provided explicit results for large current deviations in the 
totally asymmetric zero-range process (with hopping rate w n = 1) and showed how 
these can also be obtained via the application of new Bethe ansatz results for the 
totally asymmetric simple exclusion process with particle-dependent hopping rates. 
Significantly, for boundary parameters giving a well-defined stationary state (i.e., no 
condensation), we found that for low currents the fluctuations are the same across all 
bonds but for higher currents they are spatially inhomogeneous. We interpreted this 
change in behaviour as resulting from the presence of "instantaneous condensates"— 
the temporary build-up of particles on some site(s). It is thus expected to be a generic 
feature for the TAZRP with other choices of w n for which the partition function Z [see 
equation (J2J)] has a finite radius of convergence. 
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We also expect to see similar behaviour for the partially asymmetric zero-range 
process (PAZRP) with the range of validity for the general result ()26|) dependent on 
the particular choice of w n . For currents larger than some critical current we again 
expect a change in character of the current fluctuations resulting from the appearance 
of "instantaneous condensates" . The PAZRP with w n = 1 can, of course, be mapped to 
the partially asymmetric exclusion process (PASEP); in principle, it should be possible 
to solve this mapped system by Bethe ansatz (at least for a small number of particles 
corresponding to small PAZRP system size, cf. [IB] for particle-independent hopping 
rates) but a determinant form is not available. Another consequence of the bidirectional 
movement of particles is that to determine fluctuations across the Zth PAZRP bond we 
would need to solve a PASEP with L + 1 particles rather than just / + 1. 

Finally, we remark that although our results can be applied to large systems by 
taking the thermodynamic limit L — > oo, this limit does not necessarily commute 
with the long-time limit t — > oo. Our analysis thus has nothing to say about current 
fluctuations in a genuinely infinite system (taking the limit L — > oo first, followed by 
t — > oo). 
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